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Abstract 

We study the mean-field approximation of Quantum Electrodynamics, 
by means of a thermodynamic limit. The QED Hamiltonian is written in 
Coulomb gauge and does not contain any normal-ordering or choice of bare 
electron/positron subspaces. Neglecting photons, we define properly this 
Hamiltonian in a finite box [—!>/ 2; L/2) 3 , with periodic boundary conditions 
and an ultraviolet cut-off A. We then study the limit of the ground state 
(i.e. the vacuum) energy and of the minimizers as L goes to infinity, in the 
Hartree-Fock approximation. 

In case with no external field, we prove that the energy per volume con¬ 
verges and obtain in the limit a translation-invariant projector describing the 
free Hartree-Fock vacuum. We also define the energy per unit volume of 
translation-invariant states and prove that the free vacuum is the unique min- 
imizer of this energy. 

In the presence of an external field, we prove that the difference between 
the minimum energy and the energy of the free vacuum converges as L goes 
to infinity. We obtain in the limit the so-called Bogoliubov-Dirac-Fock func¬ 
tional. The Hartree-Fock (polarized) vacuum is a Hilbert-Schmidt perturba¬ 
tion of the free vacuum and it minimizes the Bogoliubov-Dirac-Fock energy. 


1 Introduction 

In Coulomb gauge and when photons are neglected, the Hamiltonian of Quantum 
Electrodynamics (QED) reads formally |HT1 ITH?1 PT/fl PHI 171 

® [ f p(x)p(y ) , , 

2 JJ l^T dxdy (1) 

where (x) is the second-quantized field operator satisfying the usual anti-commutation 
relations, and p(x) is the density operator 

, ^ _ <*[**(*), (a:)] Et=i {**(*)** (a)* - (aO^a:)*} 

P[X) 2 2 ’ [Z) 

where cr is the spin variable. In (|T|). D° = —ia ■ V + mo/3 is the usual free Dirac 
operator m (mo > 0 is the bare mass of the electron), a is the bare Sommerfeld 
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fine structure constant and p is the external potential. We have chosen a system of 
units such that K = c = 1. 

In QED, one main issue is the minimization of the Hamiltonian O- For the 
study of systems like atoms or molecules, one usually uses the Born-Oppenheimer 
approximation and treats the nuclei as fixed and time-independent external sources, 
justifying the introduction of the external field p. A global minimizer of HT p would 
then be interpreted as the free vacuum if p = 0, or the polarized vacuum when 
p ^ 0. On the other hand, a minimizer in the TVth charge sector (in a sense to be 
made more precise) would describe a bound state with particles. 

However, none of these minimization problems make sense a priori , since HF 
cannot be bounded from below. This old problem IToI fTB ; is due to the fact that 
the free Dirac operator D° has a negative essential spectrum: a(D°) = (—oo; —mo]U 
[too;oo). Our goal in this paper will be to use a thermodynamic limit in order to 
give a precise mathematical meaning to the minimization of EF in the mean-field 
approximation, that is to say when it is restricted to Hartree-Fock (HF) states. We 
shall in particular be able to define properly the Hartree-Fock global minimizer (i.e. 
the vacuum), with or without an external field <p. 

We have neglected photons in 0 - This approximation is physically justified 
when studying the free vacuum (i.e. when p = 0). In full QED, the expectation 
value of the photon field (A(x)) vanishes for the (non Hartree-Fock) free vacuum, 
see e.g. mmm- But photons should be included for a complete treatment of 
the external field case. 

We emphasize that 0 does not contain any normal-ordering or notion of (bare) 
electrons and positrons: 'F(j:) can annihilate electrons of negative kinetic energy. 
Indeed, the distinction between electrons and positrons should be a result of the 
theory and not an input. The commutator used in the formula 0 of p{x) is a kind 
of renormalization, independent of any reference. It is due to Heisenberg m (see 
also jl3l Eq. (96)] and 03 Eq. (38)]), and has been widely used by Schwinger (see 
jl%l Eq. (1.14)], 0§ Eq. (1.69)] and of) Eq. (2.3)]) as a necessity for a covariant 
formulation of QED. More precisely, the Hamiltonian HF possesses the interesting 
property of being invariant under charge conjugation since the following relations 
hold formally 

e 6’p(x) < te~ 1 = -p(x), = HT*’, 

where is the charge conjugation operator acting on the Fock space (details will 
be given later on). Notice that the use of a commutator in the same way as 0 for 
the kinetic energy term would have no effect since the Dirac matrices are trace-less: 

i J tr[fH*(ar), D°^(a;)] dx = j 'i’*(x)D 0 'l!(x) dx. 


In the mean-field approximation, one restricts the QED Hamiltonian 10 ) to 
Hartree-Fock states in the underlying Fock space. These states |f2) have the prop¬ 
erty of being totally described, up to a phase, by the two-point function P(x,y) = 
(fI|'F*(:r)\I/(y)|f2), also called the density matrix of |fl). In the usual Hartree-Fock 
theory, the energy only depends on the orthogonal projector P. This is also true 
in QED but, due to the specific definition of the density operator p(x) in 0, the 
energy will be more easily expressed in terms of the renormalized density matrix 


7 := 


P- P ± 
2 


P-1/2 


( 3 ) 


where we have introduced P^ = 1 — P, or equivalently 


7 (x,y) 


(n| 


|H). 
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The operator 7 and the associated projector P = 7 + 1/2 will be our main objects 
of concern in this paper. Notice that 7 is closely related to the Feynman propagator 
Sf taken at equal times |[7j, a widely used object in QED (for the true non-HF QED 
vacuum): 

lix, y ) = —iS F (x, y; t x = t y )/3. 

In order to give a precise mathematical meaning to the minimization of HT p in the 
Hartree-Fock class, we proceed as follows: we first define properly HF in a finite box 
in space Cl = [—L/2; L/ 2 ) 3 and with an ultraviolet cut-off A in the Fourier domain. 
For the sake of simplicity, we also use periodic boundary conditions on Cl , he. we 
work on the torus := K 3 /(LZ 3 ). Since the problem becomes finite dimensional, 
the minimization of this well-defined Hamiltonian restricted to Hartree-Fock states 
makes sense. Let us denote by El( p) the minimum energy in the presence of the 
external field ip. The main goal of this paper will be to study the thermodynamic 
limit L —> 00 , i.e. the behavior of both the minimizers and the minimum energy 
El(p) when the size of the box grows. 

Our states will always be represented by their density matrix P or their renor¬ 
malized density matrix 7 = P — 1/2, defined in @. We shall prove that, as L goes 
to infinity, the sequence of global minimizers indeed converges to a state defined on 
the whole space, which will be interpreted as the Hartree-Fock global minimizer of 
the QED Hamiltonian. 

Moreover, our method will also allow us to define the energy of these states. We 
have to study separately the free case ip = 0 and the external field case ip ^ 0. In 
the free case, we shall define properly the energy per unit volume of translation- 
invariant states; the HF free vacuum will be the unique minimizer of this energy. In 
the external field case, we will obtain at the limit the Bogoliubov-Dirac-Fock (BDF) 
energy m which is the energy of the state, measured with respect to the energy 
of the free vacuum. The HF polarized vacuum will be a minimizer of this energy. 
Therefore, the BDF theory is the appropriate model for the study of Hartree-Fock 
states in QED in the presence of an external potential. 

These limits are studied for a fixed ultraviolet cut-off A. Our objects will be 
logarithmically divergent in A. The elimination of this last divergence would require 
a renormalization , which we do not address in details here. 

Let us now describe the results of this paper in more details. 

We start with the free case ip = 0. First, we prove the existence, for each L, of 
a unique global minimizer of El(0), 7 3 = V £ ~ 1/2, which is moreover invariant 
under translations. We then study the energy per unit volume and prove that 


lim 

L —>00 


El( 0) 
L 3 


min 

7=/(—*V), 

— !/2<7<l/2 


n/)- 


(4) 


The energy T is a natural functional defined for translation-invariant operators only, 
7 = /(—iV), by 


n/) = 


(27T) 3 


/B(O.A) 


tr C 4 [D°(p)f(p)]dp- 


(2t r) 5 


B(0,A) 2 


f- r c 4 [f(p)f(q)] 

\p~q\ 2 


dpdq, (5) 


where we recall that A is the ultraviolet cut-off. The functional T indeed can be 
interpreted as the energy per unit volume of translation-invariant states. 

We prove that, when L —> 00 , the sequence ( 7 °) converges in some sense to 
the unique translation-invariant minimizer of the r.h.s. of 0 , denoted by 7 0 = 
V°_ — 1/2 where V°_ is an orthogonal projector. This state is interpreted as the HF 
free vacuum. Writing the Euler-Lagrange equation satisfied by 7 0 , we obtain the 
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self-consistent equation 


7° = 


V u = D — a 


sgn( 2 ?°) 

2 ’ 

o . 7° ( x >y) 
\x-y\ ' 


Written in terms of V°_ = 7 ° + 1/2, this equation reads 


P-=X(-oo;0)CD°), 


( 6 ) 


(7) 


which corresponds to the usual Dirac picture that the density matrix of the free 
vacuum should be the projector associated with the negative part of the spectrum 
of a translation-invariant Dirac operator. 

Let us emphasize that in many works, the free vacuum is assumed to be repre¬ 
sented by the free Dirac projector 

P- '■= X(-oo;0)( D °), 

following thereby ideas of Dirac mm (see, e.g. mmmEmEziiaEmEsi)- The true 
Hartree-Fock vacuum V_ obtained in this paper is different from P _, except when 
a = 0. The interpretation is that, contrary to P_, the self-consistent interactions 
between the virtual particles of the vacuum are taken into account via the term 
a 7 °(a; — y)\x — y| _1 in @. 

The self-consistent operator T>° defined in © has an interesting special form. 
Using the usual notation p = — *V, it can be written 

V°(p) = a ■ p 91 ^ + Pg 0 (\p\) 

\P\ 

which means that, due to ©, 


7 °(P) 


_ gi(bl) _ n n _ go(bl) q 

2blv/gi(bl) 2 + go(bl) 2 WgM 2 + goi\p\) 2 


( 8 ) 


In QED, the Feynman propagator Sf is often expressed using the Kallen-Lehmann 
representation ©3 ©HP, based on relativistic invariances. Although our model is 
not fully relativistically invariant (we discard photons and use an ultraviolet cut-off 
A) and is only defined in the mean-field approximation, our solution © has exactly 
the form which may be derived from the Kallen-Lehmann representation for the 
equal time propagator. In four-dimensional full QED, a self-consistent equation 
similar to © is well-known and used. These so-called Schwinger-Dyson equations 
□303 have been approximately solved for the free vacuum case first by Landau et 
al. in P0, and then by many authors (see, e.g., 1311230) 

In UU, Lieb and Siedentop arrived at the same equation © with totally dif¬ 
ferent arguments from ours. In particular, they did not derive the equation © as 
the Euler-Lagrange equation associated with a minimization problem, but rather 
looked for a self-consistent normal-ordering in a free Hamiltonian. They proved 
the existence of a solution of © by means of a fixed-point method, valid under a 
restrictive condition of the form a log A < C where A is the ultraviolet cut-off. Our 
proof is completely different since it proceeds by minimizing the functional T. This 
enables us to prove the existence of a solution to © without any constraint linking 
a and A. 

Let us remark that since 7 0 is translation-invariant, the associated density of 
charge , formally defined by p 1 o{x) = tr^ 7 °(a;, a;), is indeed a constant. A conse¬ 
quence of the special form © of 7 0 is that o(ir) = 0, the Dirac matrices being 
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trace-less. Therefore, in this formalism and thanks to 0. the free vacuum has no 
local density of charge, which is physically comforting. 

We know that 7 0 is a minimizer of the energy T among other translation- 
invariant operators. Since however 7 ° is for any L a global minimizer, we will 
also prove (in a sense to be made more precise later) that the limiting 7 0 is also 
a minimizer of H° under local Hilbert-Schmidt perturbations. Our study therefore 
shows that although one cannot give a meaning to the energy of the free HF vac¬ 
uum, one can give a sense of being a minimizer of the Hartree-Fock energy, either 
among translation-invariant operators, or under Hilbert-Schmidt perturbations. 

We then study the external field case ip ^ 0. We prove that the energy measured 
with respect to the free energy has a limit: 

lim {El(p) — El( 0 )} = min £bdf( 7 ~ 7 °)- (9) 

l^oo 7 _ T ° 6 e 2 , 

*- l / 2 < 7<!/2 


Here £g DF is the so-called Bogoliubov-Dirac-Fock (BDF) energy i)2}il TiA 25. ?(i ; . 
associated with the free vacuum 7 0 . It measures the energy of a HF state Q) with 
renormalized density matrix 7 , relatively to the energy of the free vacuum |flo), 
whose renormalized density matrix 7 0 has been defined previously, i.e. 

£bdf( 7-7°) = “(n|HHfi> - (n 0 |H°|n 0 )". (io) 

The BDF energy reads formally 


££df(Q) = tr[V°(V 0 + QV° + +V°_QV°_)} - j <p(x)p Q (x)dx 


a ff P Q (x)p Q (y) a 

2 JJ \x-y\ dXdy 2 


\Q{x 7 y )\ 2 

\ x - y\ 


dx dy , 


(11) 


where P q{x) = tiQi(Q(x,x)) is the charge density, a well-defined object when Q £ 
62 , thanks to the ultraviolet cut-off A E3 

For any L , we prove the existence of a Hartree-Fock global minimizer for 
7 i =Vl — 1/2. Then we show that Jl — 7° = 'Pl — 'Pl converges, in some sense, 
to a global minimizer Q = 7 — 7 0 of the BDF energy CD- Hence, the solution 
7 = 7 ° + Q of the r.h.s. of m is a Hilbert-Schmidt perturbation of the free 7 0 and 
solves the self-consistent equation 


sgn(£>) 

2 


( 12 ) 


with 


V 


x>° 

D° 


+ ap( 7 _ 7 o) * 


_ ( 7 - 7 °)(x,y) 

•I “ \x-y\ 


ip + ap 7 * 1 -r 


7 (x,y) 

\ x ~y r 


(13) 


where we have used the definition of T>° in © and p 7 0 = 0. Written in terms of 
the projector V- =7 + 1/2, (tT2ll can be written 


P- =X(-oo;0)P)> 


(14) 


which one more time corresponds to Dirac’s interpretation that the density matrix 
of the vacuum should be the projector associated with the negative spectrum of 
an effective Dirac operator. We emphasize that, due to the self-consistent 
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equation d does not depend on the reference 7 0 used for the definition of the 
BDF energy. 

Self-consistent equations for relativistic Hartree-Fock states are well-known in 
QED. We refer the reader for instance to uni Eq. (4)] which is exactly equivalent 
to Ijl4[l . and to [23 E3 [23 El ED] for related studies. In [231 EH, similar equations 
are obtained in the relativistic density functional theory. 

Notice that the BDF energy m can also be used in the free case p = 0. 
Similarly to PI I5| 1551. we shall prove in this case that its minimizer is 0 , which 
precisely means that the free vacuum 7 0 is a minimizer of H° under Hilbert-Schmidt 
perturbations. 

The Bogoliubov-Dirac-Fock energy has been introduced first by Chaix-Iracane 
in mm- They also study a mean-field approximation of QED, but start with a 
different Hamiltonian: taking the free Dirac projector P® as a definition of the free 
vacuum, they used a Hamiltonian normal-ordered with respect to P_ . In i 1 -'ll , Chaix, 
Iracane and Lions proved that their free vacuum P£ is stable when no external field 
is present. Our BDF energy m is similar to the one of Chaix-Iracane, but P_ and 
D° have been replaced by V°_ and T>°. 

The BDF model of Chaix-Iracane has been studied by Bach, Barbaroux, Helffer 
and Siedentop in [3, in which a mathematical setting without any ultraviolet cut¬ 
off is provided when ip = 0, to prove the stability of the free vacuum P°_ under 
trace-class perturbations. The case of an external field is studied as well, but the 
vacuum polarization is neglected, leading to a totally different model. 

A rigorous framework for the study of the external field case has been recently 
provided by Hainzl, Lewin and Sere in Esmg. There, an ultraviolet cut-off is 
introduced and the energy is defined on a set of Hilbert-Schmidt operators which 
are not necessarily trace-class. Starting from the Hamiltonian of Chaix and Iracane, 
it is proved that the corresponding BDF energy is bounded from below and that 
it possesses a minimizer. It satisfies an equation similar to m, but which still 
explicitly depends on the chosen reference P£ [23 Eq. ( 6 )]. We shall rely heavily 
on [23123. In particular, we shall generalize the results of m to the new BDF 
energy GU and obtain the existence of a minimizer satisfying our equation o. 
For a time-dependent study of the BDF model of Chaix-Iracane, we refer to ESI- 

Summarizing our results, we have been able to give a meaning to the Hartree- 
Fock approximation in no-photon QED, by means of a thermodynamic limit. The 
free vacuum has a renormalized density matrix 7 0 = V°_ — 1 / 2 , which is a translation- 
invariant operator, solution of ©• It is a minimizer among other translation invari¬ 
ant operators of the energy per unit volume T defined in m and as well a minimizer 
of H° under Hilbert-Schmidt perturbations. In the presence of an external poten¬ 
tial, the HF polarized vacuum is a Hilbert-Schmidt perturbation 7 = 7 0 + Q of the 
free 7 0 . The operator Q minimizes the BDF functional m and 7 solves (PI). 

This work shows that the BDF functional m is the appropriate tool for the 
study of Hartree-Fock states in QED, in the presence of an external potential. 
Although the present paper only addresses the case of the vacuum, our method 
should be applicable to the case of atoms and molecules. In that case, one would 
have to minimize the BDF energy in a sector of fixed charge N. 

As mentioned above, the theory is still divergent in the cut-off A. This phe¬ 
nomenon is also encountered in full QED, where the concepts of mass and charge 
renormalization EEI0 are used to get finite physical quantities. Mass and charge 
renormalization will usually not remove the divergences in the propagator Sf or in 
the equal time propagator 7 . This is not surprising since 7 as defined in 0 does not 
represent a physical quantity, it is not the propagator of a physical particle. The 
physical electron (at rest) is represented by the lowest energy state in the charge one 
sector (the positron corresponds to the charge -1 sector). In particular, the rest mass 
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of the physical particle is the lowest energy in the charge one sector. The physical 
electron is not just the bare particle, but includes also a cloud of virtual particles. 
The relation between the bare propagator Sf and the propagator for the physical 
particle is in the physics literature often expressed by what is called a wavefunction 
renormalization. We will discuss the concept of renormalization in Section 12.51 As 
we will explain, wavefunction renormalization is conceptually different from charge 
and mass renormalization. A completely different notion of wavefunction renormal¬ 
ization, which also leads to a different notion of physical mass, was introduced in 

EH- 

Of course, this study only applies to the Hartree-Fock approximation, the general 
case being much more difficult in our point of view. In particular, we do not know if 
the use of a normal-ordered Hamiltonian with respect to the free electron/positron 
spaces defined by V'f as proposed in eh , is physically relevant. By definition, 
the normal-ordering procedure takes a projector (i.e. a Hartree-Fock state) as a 
reference, whereas the true vacuum is known to be a non Hartree-Fock state. The 
study of a thermodynamic limit for the full QED model could be a better approach 
in the quest of a non-perturbative formulation of QED. 

The paper is organized as follows. In the next section, we define properly the 
models and state our main theorems. We start by giving a meaning to the QED 
Hamiltonian in the box Cl and formulate the thermodynamic problem mathemat¬ 
ically. Then, we define the models on the whole space that will be obtained as 
thermodynamic limits: the energy T of translation-invariant projectors and its 
minimizer 7 0 = V°_ — 1/2, and the associated Bogoliubov-Dirac-Fock model. Fi¬ 
nally, we state our main results concerning the thermodynamic limit. For the sake 
of clarity, we have brought all the proofs together at the end of the paper, in Section 
m Section^] is devoted to the discussion of the renormalization problem. 

Acknowledgment. The authors are thankful to Jan Derezinski for useful ad¬ 
vice concerning the choice of the Hamiltonian. They also thank Isabelle Catto, 
Bergfinnur Durhuus and Eric Sere for useful discussions, and acknowledge support 
from the European Union’s IHP network Analysis & Quantum HPRN-CT-2002- 
00277. 

2 Model and main results 

In this section, we present the models and state our results. The proofs are given 
separately, in the next section. 

2.1 QED energy on the torus in the Hartree-Fock approxi¬ 
mation 

Let us start by defining the no-photon QED Hamiltonian in the cube Cl = [-L/2; L/2) 3 . 
As mentioned in the introduction, we add periodic boundary conditions and there¬ 
fore work in the torus T l '■= M 3 /(LX 3 ). 

Notations 

A function if in L 2 (Tl,C 4 ) = Lp er (Ci,C 4 ) can be written as 

1>(x) = E & k y kX = (t) 1 ^ me k {x), 

fc£27i-Z 3 /L ' ' k<i_2n1 3 IL 
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where eu{x) = e ' lk ' x /L 3 / 2 and ^{k) € C 4 . We now add a cut-off in Fourier space 
and define the following finite-dimensional subspace of L 2 (Tl,C 4 ) 

•£>a := span {e^Eo. | k £ T^, a £ {1,4}} , := 2ttZ 3 /L n B(0, A), (15) 

(e (T ) cr _ li 4 being the canonical orthonormal basis of C 4 . Note that any operator 7 
acting on fj)\ has a kernel of the form 

l(x,y)= Y 7(M) e k{x)ei(y) 
k,ie rj 

where 7 (fc, Z) is a C 4 x C 4 matrix such that 7 [£ CT efc] = 7G> k) a ' }a eiE a '. 

Its density is defined as 

P-t{x) = toc*{ 7 {x,x)) = L ~ 3 Y tr C 4(7 (k,l)) e lx ' {k ~ l) . 

A translation-invariant operator T acting on satisfies T(k,l ) = g(k)6u where 
g(k) is, for any fceri, a C 4 x C 4 matrix. Denoting now 

g{x) = (2 tt/L) 3/2 ^ g(k)e h (x), 

kdTl 

we easily see that 

T(x, y ) = ( 27 r) _ 3 / 2 3 (a; - y). 

The density of a translation-invariant operator is a constant: 

Pt(z) := tr C 4 T(x, x) = (2tt) _3/2 tr C 4 p(0) = tr C 4 g{k). 


The identity on is l]( whose kernel is 


l\(x,y) = jjY Ia 


0 ik\x—y) 


(16) 


fcerT 


and whose density is Pjl = 4|T^|/L 3 . In the following, we shall denote 


r Plk 2 l r Al 
PA 2 L 3 


(17) 


Dirac operator, Coulomb potential and external field on Tl 

Recall that the free Dirac operator is defined as 

D° = —ia • V + too/3 

where mo > 0 is the bare mass of the electron and a = ( 01 , 02 , 03 ), (3 are the 
Dirac matrices m- On the torus, we use the same notation and simply define D° 
as the multiplication operator in the Fourier domain, by (a ■ k + mo/3)^, g ri. The 
same abuse of notation will be done for the operator —i V acting as a multiplication 
operator in the Fourier domain by (fc) fcer r. 

We shall use a Coulomb potential similar to the choices made in mnuEHini. 
We define Wl as being the unique solution of 

j -AW l = 4tt ( Y 6 l* 
min Wl = 0 . 





This means that there exists a constant p > 0 such that 


W L (x) = 


L 3 


E 

fce( 27 r)Z 3 /L 

\ fe# 0 


4n 

w 


ik-x 


+ pL 2 


Notice that in H3 El El ED], p is replaced by 0 (i.e. the integral of the Coulomb 
potential is assumed to vanish). Indeed most of our results are valid if p is replaced 
by any non-negative real number. However, the non-negativity of the Coulomb 
potential, which seems physically relevant to us, is usually needed in the study of 
Hartree-Fock minimizers. Moreover, our choice better mimics the behavior of l/|fc| 2 
at 0. 

In the following, we shall also consider an external potential of the form 

1 

if = an * i—y 


in the whole space, created by external sources for instance like nuclei. We will 
assume that n £ C where C is the so-called Coulomb space m 

<M) 

We do not necessarily assume in this paper that n is an L 1 non-negative function 
but the reader should think of f R3 n = Z as being the total number of protons in 
the nuclei. 

Remark 1. Notice that, since we will add an ultraviolet cut-off, the definition <GJ) 
also contains the usual Coulomb potential since the (regularized) Dirac measure 5\ 
defined by 

= (2tt)- 3 / 2 1 b(0iA) (20) 

belongs to the space C. 

For the sake of simplicity, we shall also assume that the restriction to the ball 
_B(0, A) of the Fourier transform n is a continuous function. This allows us to define 
the external potential on the torus by 

ul( x) = £ n(k)e ik *, (21) 

fe6l A 

In the whole space, the Dirac operator with external potential will be denoted by 
D v = D° — ip = D° — an * -pj, whereas we use the notation 

Df = D Q -ip L , <p L (x) := an L * W L {x) = a f n L (y)W L (x - y)dy (22) 

Jc L 

for the corresponding operator acting on ? i A . 


No-photon QED Hamiltonian on T A 

We are now able to define and compute the no-photon QED energy in T A . The 
Fock space associated with is 

( m 

A 

i= 1 


9 






(this is indeed a finite dimensional space). On J 7 ^, the creation operator HU is 
defined as usually, for k £ T{( and {1,2,3,4}, by 

A • • • A ip p ) = (e k e a ) A tpi A • • • A 4’ p . 

It satisfies the anti-commutation relation 

(23) 

In the following, we shall use the common notation ^{x) a = J2 keT L e k {x)^f k ,er for 
the second-quantized field operator. We shall also need to use, for some f £ Sjj^, 
the notation 

/«\ 3/2 

E f( k )^k,a (24) 

' ' fcer^ 


when / = (^) 3/2 Efcerj fM)e k . 

The no-photon QED Hamiltonian is the well-defined operator acting on the 
finite-dimensional space ■/-]{', 


= f '^*{x)D a ^{x)dx — f tfiL(x)p(x)dx 

Ml M l 

+ -z[[ p{x)p(y)W L (x-y)dxdy (25) 
1 J J(t l ) 2 


where a is the bare fine structure constant and 
P(x) = ^tr[**,$](a:) 

= ^E (**(*)**(*)*-*(*)**>)') 

(7 

= ^E E ( 26 ) 

a MerJ 

The charge-conjugation operator *€ is the uniquely defined (up to a phase) unitary 
operator acting on the Fock space J 7 ^ (see, e.g. (T7I Prop. 2.1]) such that, for any 
/ e Si 

wtjye - 1 = ®(c/)* and wMytf - 1 = ^(c/), 

where C is the charge conjugation operator acting on defined by Cf := ifdazf, 
and H >{f) has been introduced in (1241) . It is then easy to see that the following 
relations hold 

tfpM)^- 1 = -p{x), = H" v . 

Using the CAR and the formula we obtain 

P(x) = V{x)*{x) - pi- (27) 
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where p\ is defined in 1171) . Inserting (EZ|), we obtain for the last term of (I251i 


7// p(x)p{y)W L (x-y)dxdy 
1 j J(Tl ) 2 

= 77 ff ^*{x)^(x)^*(y)^{y)W L (x-y)dxdy 

z J J(t l ) 2 

~ a P\ ff ^*{x)^(x)W L {x - y)dxdy + ff W L (x-y)dxdy 

JJ(Jl ) 2 1 l ) 2 

= a V W L (x i -x j ) + ^ ff I^{x 1 y)'H*{x)^{y)W L {x~y)dxdy 


l<i<j 


a PA ff {x)^[ x)Wl{x — y)cfa ;dy + 

J J(T l ) 2 


<P L a? 


// W L (x-y)dxdy, 

Ti) 2 

(28) 


where we recall that I^{x } y) is defined in IItill . Therefore, computing the energy 
of a Hartree-Fock state with density matrix P(x,y) (an orthonormal projector on 
,£))(), we obtain from CP and (11281) 


£i™( p ) = tr(£>°P) - f <p L (x)(p P (x) - p\) dx 

+ t // Pp{x)pp(y)W L (x-y)dxdy-^- \P(x,y)\ 2 W L {x - y)dx dy 

z JJ(Tl) 2 z JJ(Tl) 2 


a 

+ - 


'(J l) 


tr c 4 (^a (x,y)P(x,y)) W L (x - y)dxdy 


a PA f f pp(x)W L (x - y)dxdy + a ^ p J^ ff W L {x-y)dxdy. (29) 

l) 2 Z JJ(Jl) 2 


Using pjL , 2 {x) = p\ and tr (D°lfc) = 0, we infer 


£?,™( p ) = £Z(P - /a/2) + | ff \lk(x,y)\ 2 W L (x - y)dxdy 


(30) 


where £)/ is the usual Dirac-Fock energy on the torus 


= tr(D^ 7 ) + ^P L (p 7 ,p 7 ) - ^ ff Y){x,y)\ 2 W L {x - y)dxdy 
z z J J(Jl) 2 


(31) 


defined for any self-adjoint operator 7 acting on and with 
D L (J,f)= ff f{x)f{y)W L {x-y)dxdy =± £ \f(k)\ 2 wL(k) > 0, 

JJ(Tl) 2 L 


ke(2nZ 3 )/L 


p-r(x) = tr C 4 j(x,x). 

The last term of CP is a constant which behaves like 

tt // \I^{x,y)\ 2 W L (x-y)dxdy ~l^oo dx (32) 

° J J{t l ) 2 z /r 3 fI 

where 5a £ L 2 (R 3 ) D L°°(R 3 ) is the Fourier inverse of (27r) _3 / 2 l B ( 0 ,A) already 
defined in pof . Because of lt32l) . this term shifts the limit of the energy per unit 
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volume by a constant and disappears when looking at differences. Hence, it will not 
play any role for the study of our thermodynamic limit, and we discard it for the 
rest of the paper. We therefore study the minimization problem 


E L (ip) := inf 


^( 7 ), 


7<£ £(£a)» 7* =7, 




1 T l 

< 7 < iA 

2 - ' “ 2 


(33) 


(C(Sjjf) denotes the space of all linear operators acting on Sjjf}. Notice that we have 
extended the set {P — /£/2, P orth. projector} to its convex hull 


<?A = 


7 e £(«a), 7* =7, 


tL tL 

< < 4 

2 “ 2 


as this is usually done in Hartree-Fock type theories im It will be proved that this 
does not change the minimizer. With this change of variable, our goal will be to 
prove 


• without external field (ip = 0): 


1. the existence of a unique minimizer = P° — 7^/2 for El(0), when L 
is large enough; 


2 . 

3. 

4. 


that the energy per unit volume 


El{ 0 ) 
L 3 


has a limit as L 


oo; 


that converge in a certain sense to a translation-invariant limit 7 0 = 
P°_ — 1/2, which will be the HF free vacuum (P°. is a projector); 

that the limiting 7 0 is the unique minimizer of the energy per unit volume , 
a functional defined only for translation-invariant operators acting on the 
whole space; 


• in the presence of an external field (<p / 0 /: 

5. the existence of a minimizer 7 ^ = Vl — /a/ 2 for £x(</?); 

6 . that the energy difference E^iip) — El{ 0) has a limit as L —> 00 ; 

7. that 7 a converges in a certain sense to an operator 7 = P_ — 1/2 which 
will be interpreted as the HF polarized vacuum (V- is a projector); 

8 . that 7 ^ 7 ° is a minimizer of the Bogoliubov-Dirac-Fock energy, which 
measures the energy with respect to the (infinite) energy of the free 
vacuum 7 0 . 


These questions are rather common in the mathematical study of thermody¬ 
namic limits mmmm- Before we answer them, we have to define properly the 
variational problems obtained in the whole space. In the next section, we define 
the energy per unit volume of translation-invariant operators acting on the whole 
space when no external field is present; its minimizer will be the free vacuum 7 . 
In Section IOI we define the associated Bogoliubov-Dirac-Fock model properly and 
prove the existence of the polarized vacuum 7 in the presence of an external field; 
this will be only an easy extension of the work already done by Hainzl-Lewin-Sere 
in [2SII2S1 Finally, we answer questions (l)-( 8 ) in Section HOI 


2.2 Definition of the free vacuum V° 

Let us now define the models in the whole space. The cut-off is implemented in the 
Fourier domain by considering the following Hilbert space 

fiA := jt/ € L 2 (R 3 ,C 4 ) I supp^ C 5(0, A)|. (34) 
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In this subsection, we consider the case where no external field is present, n = 0. 
We want to define the energy per unit volume of a translation-invariant operator 7 
acting on f)A and such that 


/a , /a 

-< 7 < —, 

2 - ' - 2 


I A being the identity on Ioa- Such an operator acts as a multiplication operator in 
the Fourier domain and can therefore be written 7 = /(—iV) or, formally 7 ( 2 , y) = 
(27r) _3 / 2 /(a: — y), with / belonging to 


Aa 



L 00 (5(0, A), 64 (C)), 



Here 54 (C) is the set of 4 x 4 self-adjoint complex matrices. Notice that its density 
of charge is a well-defined constant: 


P 7 {x) = tr C 4 j(x, x) = (2tt) 3/2 tr C 4 /( 0) = (2tt) 3 f tr c *(f(p)) dp. 

Jb( o,a) 

The energy per unit volume of such a translation-invariant operator 7 = /(—zV) is 
then defined in terms of / £ Aa by 


nn 


7 ^x 3 / tr C 4 [D°(p)f(p)}dp 
(^ 7r ) Jb(o,a) 


(2tt) 5 IIb(o,a) 2 


trc *[/(p)/(g)] 

b-# 


dpdq. 


(35) 

To give a hint why we consider the energy T, let us fix some operator 7 = /(—zV) 
where / £ .4 a- To simplify the discussion, we also assume that trc*(/(&)) = 0 for 
any fc £ 5(0, A). To avoid any confusion, we denote by 7 L the operator acting on 
.£)a associated with 7 , and which is just the multiplication operator by (/(fc))fcer L 
in the Fourier domain. Its energy is given by m and can be expressed as 


£°( 1 L )= tr C 4(T>°(fc)/(fc)) -I ' 3 2 ® 

fcer£ n) 



\f L {x)\ 2 W L (x)dx 


where /l(x) := (2n/L) 3 ^ 2 f(k)e k (x). Notice that we have used 


PjL = L 3 ^2 tr C 4 /(ft) = 0 

fcer A 


since by assumption trc 4 /(ft) = 0 for any k £ 5(0, A). It is easy to see that 

£° L h L ) 


lim 


= T(f)- 


L—>00 L 3 

It is therefore natural to define the free vacuum as a minimizer of T and we introduce 


E r := inf {T(/) | / £ ^ A } ■ 


Theorem 1 (Definition of the free vacuum). Assume that 0 < a < 4/77 A > 0 
and that uiq > 0. Then T possesses a unique global minimizer f on Aa and any 
minimizing sequence ( f n ) C .4 a ofT converges strongly to f in L 2 (B(0, A), 54 (C)). 

The associated renormalized density matrix 7 0 := /(—zV) is a translation- 
invariant operator satisfying the self-consistent equation 


0 sgn(5°) 

7 =-o-■ 


V u = D° - a 


7 °{x,y) 

\x - y | 


(36) 
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or, written in terms of the translation-invariant projector Vf = 7 0 + /a/2, 

V°- = X( —oo;0) (x>°) . (37) 

Moreover, T>° takes the special form, in the Fourier domain, 

V°(p) = a-uj p g 1 (\p\) + go(\p\)P (38) 

where ui p =p/\p\ and go, g\ £ L°°([0; A],R) are such that 

gi(x) > x and go{x) > mo, (39) 

moSiW < 9 o{x)x (40) 


for any x £ [0; A), and therefore 

ml + \p\ 2 < \V°(p)\ 2 < 90 ^\ ml + \p\ 2 ). (41) 

mo 

Finally , V°(p) £ 0,A)) C C“(B(6jV)). 

The proof of Theorem |T| is given in section IQ 

This result is a generalization of a work by Lieb and Siedentop. In eh, the 
equation 111 (il l is solved by means of a fixed point method, only valid under a con¬ 
dition of the form a log A < C. Thanks to the variational interpretation using the 
function T, we have been able to prove the existence of a solution of without 
any constraint linking a and A, and by means of a completely different proof. Our 
solution coincides with EH when a log A < C, and the properties of V° stated in 
Theorem Q are exactly the ones which have been proved in EH- The interpretation 
of lllfill given in EH however does not seem to be the same as ours. 

We notice that a self-consistent equation similar to lllfill (written in terms of 
the four-dimensional Green function of the electron) has been approximately solved 
first by Landau et al. in DEI, and then by many authors (see, e.g., EH 1231EI)- 
They use an ansatz analogous to foBt (see, e.g., [T| Eq. (1)] and Eq. (4.1)]). 

Remark that the density of the vacuum 7 0 = /(—zV) also vanishes, since lllfill 
and ISHt mean 


Kp) 


gi(bl) ,_ go(bl) q 

2\/«7i(|p|) 2 + «7o(bl) 2 ' " W9i(\p\) 2 +9o(\p\) 2 


and therefore 

p 7 o = (2?r)“ 3 [ tr C 4 f(p) dp = 0, 

7.8(0,A) 

the Dirac matrices on and /? being trace-less. 

Remark 2. Adapting the proof given in Section I, 1 ?.M it can be shown that when 
mo = 0, the functional T also has a unique minimizer f(p) = —V°(p)/(2\V°(p)\), 
with T>°(p) = a ■ u> p gi(\p\), gi being given by \j /[ Eq. (26)/. 

Remark 3. The energy T satisfies a nice scaling invariance property (already used 
in mb Namely one has, with an obvious notation, 

T mo , a Af) = A- 4 T Am0iQ , AA (/(-/A)) 

which implies that 

E T (m 0 ,a,A) = A“ 4 £ T (Am 0 ,a, AA) and (V°_)xm 0 ,a,XA(p) = (V°_)m 0 ,a,\(jpl X)- 
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2.3 The Bogoliubov-Dirac-Fock theory based on V°_ 

In the previous section, we have defined the HF free vacuum, whose renormalized 
density matrix is 7 0 = V°_ — 1 a/2. We now assume that n / 0 (recall that p = 
an* 1/1 • |) and define the Bogoliubov-Dirac-Fock model mini ini 0123113, based 
on this new reference. Our goal will be to show that the results by Hainzl-Lewin- 
Sere Gam can be extended to this case. We refer the reader to G3 EU123 123 
for a detailed presentation of the BDF model. 

The BDF energy reads |25l : 25j 


^bdf (Q) '■= tr P° ( D °Q ) + ^ D (PQ, PQ ) ~ a D{pQ, n) - | jdx dy 

(42) 

where 

Q e Qa := {q £ ©f° (Sa), -V°_ <Q< V° + , PQ £ C j, (43) 

c = {f | D(f,f) <oo} 

and _ 

'P- 

The space ©W^a) is introduced in jzsj Sec. 2.1]. It contains all the Hilbert- 
Schmidt operators Q € © 2 (hA) which are such that Q ++ = V\QV\ and Q~~ = 
'P^Q'P'f are trace-class (£ ©i(bA))- The P°-trace of Q is then defined by 

tr-po Q := tr(Q ++ + Q ). 

Finally, P q{x) = trc 4 Q(x, x) is the charge density associated with Q , a well-defined 
object in Fourier space thanks to the ultraviolet cut-off 123: 

PQ{k) = (2tt)“ 3/2 [ tr C 4 (Q(p+ k/2,p- k/2)) dp. 

Jb( o,a) v 1 

Notice that compared to [13 H3 CH El E3 ESI , we have not only replaced by 
V°_ , but also D° by V° in the definition 14211 of the BDF energy. 

The interpretation of gll) is that £b DF (7 — 7 °) is the energy of the HF state 
7 , measured relatively to the energy of the free vacuum 7 0 . This statement will 
be made more precise in the next section, in which we prove that — E^{ 0) 

converges to the minimum of the BDF energy. Notice when 7 satisfies —Ia/2 < 
7 < /a/ 2 , then one has — V/_ < 7 — 7 0 < V+, justifying the constraint imposed on 
Q in the definition of Qa in ©. 

We may now state a result similar to G3 Theorem 1]: 

Theorem 2 (The BDF energy is bounded-below). Assume that 0 < a < 4/-7T, 
A > 0, mo > 0 and that n £ C. 

1. One has 

VQ £ Qa, i ^bdf( < 3) + — D(n, n) > 0 
and therefore £ is hounded from below on Qa- 

2. If moreover n = 0, then £g DF is non-negative on Qa, 0 being its unique 
minimizer. 
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Proof. The proof is the same as in D3HS1ES], using Kato’s inequality 


due to TO . 


JL < 1\D°\ < -\V°\, 

x ~ 2 1 1 “ 2 1 


□ 


The interpretation of the second part of TheoremElis that 7 0 is not only a mini- 
mizer among translation-invariant projectors fSection l 2 . 2 l l. but also among Hilbert- 
Schmidt perturbations. 

Let us now define the BDF ground state energy in the presence of an external 
field: 

L'BdfCvO = inf^BDF- (44) 

Wa 

The existence of a minimizer is obtained by a result analogous to esi Theorem 1]: 

Theorem 3 (Definition of the polarized vacuum). Assume that 0 < a < 4 / 7 T, 
A > 0, too > 0 and that n £ C. Then £b DF possesses a minimizer Q on Q\ such 
that 7 = 0 + 7 ° satisfies the self-consistent equation 


7 = 


sgn(£>) 

2 ’ 


1 7 (x,y) 


V := D { + a (piy — n) * 7—7 — a . 

M \x-y I 

or, written in terms of the projector P - := 7 + /a/2 = Q + Vf, 

'P- = X(-oo;0) (P) ■ 

Additionally, if a and n satisfy 


0 < a— { 1 — a 


f 1 ~ + ’ 1/S2ll/ “) Wc } 


(45) 


(46) 


(47) 


then this global minimizer Q is unique and the associated polarized vacuum is neu¬ 
tral: tr-po ( Q ) = 0. 

Proof. The proof works analogously to Hi Proof of Theorem 1]. The only modifi¬ 
cation is ESI Lemma 1]: one proves that || |2?°|] || e = 0(1/7?) using the 

regularity of T> {> (p) stated in Theorem |T| □ 


2.4 Thermodynamic limits 

Let us first state a result for the thermodynamic limit when n = 0, and which 
answers to the questions (l)-(4) of Section TOl We recall that E r , T, /, 7 0 and 
V°_ are defined in Section E3 

Theorem 4 (Thermodynamic limit with no external field). Assume that 
0 <a< 4 / 7 r,A> 0 , and that mo > 0. Then for L large enough, £° possesses a 
unique minimizer 7 ° = — on Q\, where T 3 / is an orthogonal projector. It is 

translation-invariant, 7 ^ = /°(—iV) (or equivalently 7 ° ( 2 , y) = (2n)~ 3 / 2 f 3 (x—y)). 
One has 

lim El ^ = E T = min{T(/) | / £ 7l A } (48) 

L —too Li 

and 

lim || 7 £ - 7° Us (TO = r lim SU P |/°( fc ) “ /( fc )| = °> (49) 

L —>00 u °° L—>00 

lim ||Vl - V°_ | s (T } = lim sup \V° L (k) - V°_(k)\ = 0. (50) 

L -*OO 11 u °°L->oo Cf-pL 
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The proof of this result is given in Section IT til 

It is interesting to see how this result would change if another choice was made 
for the definition of the no-photon QED Hamiltonian on the torus, like 


^ := [ A>*(x)D°At(x) dx + £ 
Jt l 2 



ty*(x)'fy(x)^*(y)'i’( u)Wl(x — y ) dxdy 


or 


Hi := + E W L ( Xi - Xj ) 

i> 1 1 <i<j 

= f ^*(x)D^fS(x) dx + ^ ft (x)^*(y)^(y)^ (x)Wl(x — y) dx dy. 

JTl 2 JJ(Tl ) 2 

In both cases, the thermodynamic limit is dramatically changed since we can prove 
the following result: 

Theorem 5 (Thermodynamic limit for other Hamiltonians). Assume that 
K is replaced by Hi or Hi in the previous study, that a > 0, A > 0, mo > 0 and 
n = 0. Let us denote by 0 < 0 a < a, self-adjoint operator acting on which 
is the one-body density matrix of a minimizer of the QED Hamiltonian, when it is 
restricted to translation-invariant Hartree-Fock states. Then 

lim p eL = lim L ~ 3 Y] tr C 4 (Q L (k)) =0. (51) 

feerj 

Therefore, writing &L(x,y) = (27r) _3 / 2 ^(a: — y), 


L li ^ 0 ^ i ^ 00 (R 3 ) = °' ( 52 ) 

In other words, when the Hamiltonian Hi or Hi is chosen instead of Hi, the 
thermodynamic limit as L —> oo is trivial. This phenomenon is due to the term 
Dl(p, p) which behaves like L 5 and not like L 3 and which plays the role of a penal¬ 
ization: the limit necessarily has a vanishing density. This shows the usefulness of 
the commutator in the definition of p(x) in (T2TT1) or In that case, one obtains 

lim p o = lim p (r o tl /2) = lim ( p v o - p L J 2) = 0 
L— s-oo L L— s-oo v L a/ / L—►oo L 

8Jld p^y o = p^p o —"pOj/2 = 0. 

We now state a result for the thermodynamic limit when n / 0; this answers 
the questions (5)-(8) of Section ITll 

Theorem 6 (Thermodynamic limit with external field). Assume that 0 < 
a < 4/-7T, A > 0, mo > 0 , n £ C and that n is continuous on 5(0, A). Then for 
any L, possesses a minimizer 7 l = Vl — Ia/ 2 on where Vl is an orthogonal 
projector, and one has 

lim {E L (ip) - El(0)} = Ebdf{t) = min {£^ DF (Q), Q € Qa} . (53) 

L —►oo 

Moreover, up to a subsequence, Ql(x,p) := (yz, — 7 ^)(x,y) = (Vl — V^)(x,y) 
converges uniformly on compact subsets of R 6 to Q(x,y), a minimizer o/£g DF on 

Qa- 

This result justifies the Bogoliubov-Dirac-Fock theory for the study of the po¬ 
larized vacuum in no-photon QED restricted to Hartree-Fock states. The proof is 
given in Section 1331 
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The study of the thermodynamic limit for atoms and molecules (i.e. when the 
minimization is restricted to a specific charge sector) is clearly out of the scope of 
this paper. Let us however indicate briefly what should be obtained in that case. 
We fix some integer N and define the minimization problems in charge sectors: 

£fM:=inf{^( 7 ), 7 e^, tr( 7 ) = JV}, (54) 


Ebuf(v>) ■■= mf{^ DF (0) | 0 

A minimizer Q = 7 — 7 ° of EU), if it exists, is 
of the form m. 

sgn (Dj - fi) 

7 “ 2 

:= D° + a (p-y — n) 

p being an Euler-Lagrange multiplier due to the charge constraint, interpreted as a 
chemical potential. Written in terms of V- = 7 + /a/2 , 15011 reads 

'P- = X( —oo; 0 ) {'Dj ~ h) = X(-oo;aO (^7) • ( 57 ) 

We refer the reader to d and |25l Remark 6] for comments in connection with the 
Dirac-Fock equations m. 

We now conjecture that 

lim {< (ip) - E L ( 0)} = Eg DF (ip) (58) 

L —>^oo 

which would also justify the Bogoliubov-Dirac-Fock approach for the minimization 
in charge sectors. 


e Qa, ti>o (Q) = AT |. (55) 

a solution of a self-consistent equation 


1 7 (x,y) 

M \x-vY 


(56) 


2.5 The Renormalization Problem 


Throughout this paper, the cut-off A has been considered as a fixed parameter. 
However, most of the objects that we have obtained, like for instance the density 
matrix 7 , are divergent in A. This phenomenon is also encountered in full QED. 

To give an example of this, let us consider <7o(0), where go has been defined in 
Theorem El Writing the self-consistent equation 63, it can be seen that Eq. 
( 20 )] 

/ n \ , « f 1 flo(M) , 

50 (0) = m 0 + -T~2 / M2 / n 

4? r J B ( 0 ,A) \l\ \/gi(\q \) 2 +5o(|5|) 

Therefore, using GD, one obtains 


5o(0) > 


mo 


a 

47T 2 


m 0 go(\q\) 


’b( o,a) M 2 5o(M)\/M 2 + m n 


mo 1 + 


47T 2 J B ( o,A) kl 2 y/\q\ 2 + ml 


m 0 ( 1 + 


a 


arcsinh(A/mo)^ 


(59) 


which shows that go{0) is divergent as A —> oo. 

In regular QED, the divergences of the (appropriately defined) physical mea¬ 
surable quantities are usually eliminated by means of a mass and a charge renor¬ 
malization. The main idea is to assume that the parameters a and mo appearing 
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in the theory are indeed bare parameters which are not physically observable. The 
physical parameters are assumed to be functions of a, toq and the cut-off A 


Oph — rTph(o, rriQ , A), m p h — r7i p h(o:, mo, A) 


and equal the physical values obtained in experiment. These functions should be 
inverted in order to express the unknown bare quantities in term of the physical 
quantities 

a = a(a p h,m p h, A), m 0 = m 0 (a ph , m ph , A). (60) 


Using these functions, one expects to remove (in some sense that needs to be pre¬ 
cised) all divergences from physically measurable quantities. 

Mass and charge renormalization however does not remove all divergences in the 
theory. Certain quantities, e.g. the bare Feynman propagator Sf (either at equal 
times or at general space time points), are still divergent. The expectation is that 
all these divergences cancel in physically measurable quantities and that they are 
therefore of no real relevance in formulating the theory. 

Although there is no real need to do this, it is often convenient to introduce 
a renormalization of the bare Feynman propagator Sp- This is referred to as a 
wavefunction renormalization. In full QED m it is claimed that the divergence in 
the Feynman propagator may be removed by a multiplicative renormalization and 
that the renormalized propagator has the same pole near mass shell in 4-momentum 
space as a free propagator corresponding to a particle with the correct physical mass. 

Note that in practice, this theoretical renormalization procedure is always used 
to justifying the dropping of the divergent terms obtained at each order of the 
perturbation theory jl ftj . For this fact to be true, it is particularly important that 
renormalization can be expressed by means of multiplicative parameters in front of 
the different propagators lia- 


In Hartree-Fock QED, it is not clear at all if the usual renormalization program 
of QED can be applied (especially when photons are not included). In m , a renor¬ 
malization of the charge is proposed but it seems to be only valid if the exchange 
term is neglected. In m p. 194- 195], it is argued that mass and charge renormal¬ 
ization is alone not enough to completely remove the divergences of the HF theory 
by means of multiplicative parameters. 


In ^j, it is proposed, among other possibilities, to write the operator of Theo¬ 
rem ^ as 


V°(p) 


■91 (M) 

\p\ 


• p + 


HffoOl) 
91(\P\) 1 J 


and to interpret < 7 o( 0 )/< 7 i( 0 ) as being the physical mass m p h and gj( 0 ) as being the 
wavefunction renormalization constant. This proposal does not at all correspond 
to the procedure described above. Indeed, multiplying T>° by a constant does not 
change the equal time Feynman propagator 7 0 , since 7 0 = — sgn(2?°)/2. 

We believe that the correct solution is rather to define the quantities a p h and 
TO p h for a free Hartree-Fock electron (strictly speaking, photons should be included). 
Namely, we define the physical mass as the minimum energy of a free electron which 
means, using (|55| . 

m ph '■= -Ebdf(O)- (61) 


Since an electron never sees its own field but interacts with the Dirac sea, we believe 
that £’g DF (0) possesses a nrinimizer Q = 7 — 7 0 , i.e. that an electron can bind alone 
when interacting with the Dirac sea (this has been first proposed by Eric Sere). 
According to m , the operator 7 would then be a solution of the self-consistent 
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equation 


7 = 


sgn(£>7 - p) 


1 


T>j := Z? u + ap y * -—- — a 


7 (x,y) 


(62) 


\x-y r 

where p is a chemical potential chosen such that tr-po (7 — 7°) = 1. Formally, we 
could then define a p h by 

\-i\ 


a p h := a lim |x| (p 7 * |■| 3 ) ( x ). 


(63) 


If py is smooth enough, for instance p y G W 1 '°°{B(Q, A)), this is equivalent to 

a ph := (27r)' 3/2 ap^(0). 

If 7—7° were trace-class, one would have (27r) _3 / 2 py (0) = tr(7—7 0 ) = f R3 py_ 7 0 = 1 
and a p h = cx. But in the present case, it is known that when 7 is a solution of the 
self-consistent equation El such that p 7 7^ 0, then 7 — 7 0 is never trace-class 
El El ES], Expanding El to first order in the bare a, this should rather lead to 
a formula of the form ra 

“ tl( “' ra °' A) = l + Sl,°g(A/m„) + ° (a) - 

the usual formula in QED El Eq. (7.18)]. 

Both El and El would define m p h and a p h as extremely complicated non¬ 
linear functions of a , mo and A. A challenging task is to study the finiteness of 
measurable quantities like for instance the energy of an electron in the presence of 
an external field E , g DF (<p), when a p h and m p h are fixed to be the observed physical 
quantities. 


3 Proofs 


3.1 Preliminaries 

We shall need the following lemma, which summarizes the properties of Wl , defined 
in itTbl) : 

Lemma 1. The potential Wl satisfies the following properties: 

1. there exists a constant C such that \)2\ 


W L {x) 


1 


L°°(C l ) 



(64) 


2. (Kato-type inequality) for any mo > 0 and L large enough, there exists a 
constant C^(mo) such that, on Sj%, 


W L < Ck(m 0 )\D°\, 


lim C A {m 0 ) = C A (m 0 ) 

L—>00 


( 65 ) 


where C A [mo) < 7t/2 is the optimal constant for the inequality 


wr < C A (m 0 )\D° 


on S) a. 
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Proof. We know from E2 Eq. (112)] that C := 
other hand, we have the scaling property 



< oo. 

L°°{Cx) 


On the 


w l(i) = l Wl (f) 


which implies that 


W L - — 


L°°(Cl) 


c 

< — . 

- L 


For a fixed bare mass mo > 0, let us define (we do not mention Too for simplicity) 

Cfc := sup f Wl(x) | ip(x) 1 2 dx < oo, 

Jc L 


V’GflA JC L 


C A ■■= sup / 
tp&s) a Jr- 


IF(F)| : 


a Jk 3 Fl 


■dx < oo. 


A proof that is finite can be found in m ■ Let us consider a function ip £ Sja 
such that ip is continuous on the ball B(0, A) and 


\ip{x)\‘ 


■dx — Ca 


<e, (\D°\ip,ip) < 1. 


Jr 3 FI 

We now approximate ip by a sequence of functions in by defining 

( 2 t r ) 3 / 2 


E V#K fc ' x - 


fcerf 


We have 


Pl\\L°°(C l ) ^ 


(27r) 3 / 2 1 r 


Al 


L 3 




L“(B(0,A)) 


IVtM 


A(27r) 3 / 2 |r^| 
A-(Cz,) S ■ 


L 3 


iP 


L°°(B( 0,A)) 


Therefore (iPl)l is bounded in W 1,0O (R 3 ) and ipL —* ip uniformly on compact 
subsets of R 3 . On the other hand, 


^\\l H C L ) ^ ' L 3 


(27r) 3 |ri 




L”(B(0,A)) 


which is also uniformly bounded in L. Using now 


f \ip L (x)\ 2 W L {x)dx - f EMLrfa; 

Jc L Jc l Fl 


^ IIFaIIl^Cl) 


W l - 


< 


L°°(C l ) 


C 

L 7 


we easily deduce that 


lim [ \ip L (x)\ 2 W L (x)dx = [ 
L ^°° Jc L JR 


\ip{x)\‘ 


-dx. 


r R3 F 


On the other hand, 


(\d°\iPl,iPl) = ('=f) E (\D°(k)\m,m) 


L I 


c 4 


fcerf 
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which converges to (\D°\ijj,ip) as L —> oo, and therefore ipi) < 1 for L 

large enough. This shows that 

lim inf Cfc > Ca + e 

L —► oo 

for any e > 0 and thus lim inf l->oo Ca > Ca- 

Since is finite dimensional, there exists for any L an optimal function /a G 
satisfying (|-D°|/l, /l) < 1 and J Cl WL(x)\fL(x)\ 2 dx = Cfc. Thanks to the cut-off in 
Fourier space, (/l)l is bounded for instance in H 1 ^^) and therefore /a —> / G $)a 
in Lf oc (R 3 ), up to a subsequence. Passing to the limit like above, we obtain 

lim [ \f L (x)\ 2 W L (x)dx = [ dx 

-tR 3 R 

and (\D°\f,f) < 1. This shows that 

lim C L a < [ dx < C A . 

L ^°° J R3 \x\ 

As a conclusion, we have 

lim Ca = Ca < 7r/2 

L—>00 

which ends the proof of Lemma ^ Cl 


3.2 Proof of Theorem ID 

The functional T can as well be written 


T(f) = 


(2tt) 3 


'B(O.A) 


tr C 4 [D°(p)f(p)]dp- 



which easily shows that T is continuous for the weak topology of L 2 (B{ 0, A), 54(C)). 
Since A a is bounded and closed in this topology, T possesses a minimizer on Aa- 
In order to show the uniqueness and the mentioned properties of this minimizer, we 
shall now minimize T on a special subset of Aa- 


• Step 1: Minimization of T on a restricted set. Let us introduce the following 
subset of Aa 


Ba := {/ G Aa | 3(/ 0 ,/i) G G, f(p) = a • w p fi(\p\) + /o(|p|)/3} (66) 


where lo p = p/\p\ and G C L°°([0; A], R) 2 is the set of pairs (/o,/i) satisfying the 
following properties 


Vp G B(0, A), 


\p\ 


'B(O.A) 


fo < 0 and fi < 0, 

/o 2 + /i <1/4 

/o(M) j / 
-2 dq < mo 






B(0,A) 


\v~qV 


dq. 


(67) 

( 68 ) 

(69) 


Notice that G is a bounded convex set, closed in the weak topology of L 2 ([0; A], R) 2 . 
Therefore, Ba is also a closed subset of Aa- As a consequence, we deduce that T 
possesses a minimizer / on 5 a • Let us now show that / satisfies the same Euler- 
Lagrange equation as global minimizers of T on Aa- 

Indeed, considering perturbations of the form (1 — t)f + tg for any g G Ba, 
t G [0; 1], one easily shows that / is a minimizer of the following energy on Ba 

Tfid) = [ tr C 4 (V°(p)g(p)) dp, (70) 

Jb( o,a) 
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i.e. 7 'f(g) > Tf(f) for any g £ Ba ■ Here T>° is the mean-field operator associated 
with / 


V\p) 


D°(p) - 


a 

2 i r 2 


m 


■dq 


a-u)p\ \p\- 


<B( o,A) b~ ?l 2 
a /' tJp ■^g/i(Jgj) dg 


27T 2 


f B(0,A) 


a 


Ip-9 I 5 


+/3 P“2^. 

a • w P 5i(|p|) +/3p 0 (|p|), 


S(0,A) 


Mill 

b- b 2 


dq . 


(71) 


with 


and 


Si(bl) = bl - 


a 

27T 2 


r B(0,A) 


• ^g/i(kl) 

\p~q \ 2 


dq 


= bl — — / riQi (l (— + rr) ) h{v)dv. 


/0 


bl \ 2 v« bl 


5o(bl) = m o - ^2 


B(0,A) 


/o(lb) 

b-b 2 


dq 


/ A A<3o(1 (M + JL 


7T 


F 0 


bl 


V 


bl 


fo(v)dz 


(72) 


(73) 


Here Qo and Qi are positive functions on (1; oo), defined in gTJ Eq. (24) and (25)], 
arising from the integration of the angle as shown in M 

Due to the assumptions of IP . we deduce that <?i(|p|) > |p| and <7o (bl) m o, 

which implies 

|Z> 0 (p)| > \J\p \ 2 +mg > m 0 > 0. 

Therefore T>°(p) is always invertible. This now easily implies that the unique global 
minimizer of Tj defined in CZ3, on A a, is 


I\P ) = 


sgn(£>°(p)) 

2 


P°(p) 

2|2?°(p)|- 


We now show that indeed f £ Ba, which will imply f = /. To this end, we just 
have to prove that 


fo(p) 


9o(p) 

2|2?°(p)| 


and f[(p) 


9i(p) 

2\V°(p)\ 


satisfy the additional properties lIHTI) and 11691 ) . For 1671 ) . this is an immediate 
consequence of the fact that Qo and Q\ are non-negative on (l;oo), whereas (Ibil is 
proved in PH . Theorem 2]. 

As a conclusion, we have found a minimizer of T on the restricted set Ba which 
is a solution to the self-consistent equation 


= sgn(£>°(p)) 


where T> 0 (p) is defined in CD- Notice that by construction, V°(p) satisfies the 
properties l!3§ll . (TTHll and iTTTTl . Of course, 7 0 := /(—zV) and V°_ := 7 0 + /a/2 are 
respectively solutions of the self-consistent equations m and E3. 
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• Step 2: the minimizer f of T on B\ is its unique global minimizer on „4a. We 
now show that we not only have solved the self-consistent equation m, but that 
/ is the unique global minimizer of T. 

To this end, we compute, for some / £ 


T(f)=T (/) + 7^T3 ( [ t r c 4 [ V°(p)Q(p)]dp - | , 

(2tt) j \Jb(o,a) 2 -« 3 1*1 / 

where Q(jp) = f(p) — f{p)- Since — Ic 4 /2 < f(jp) < Ic 4 /2 due to EHt and f(p ) = 
(P'ffp) — V+(p))/ 2, we see that 

Vp g -8(0, A), -V°_(p) < Q(p) < P°(p). (74) 

We now adapt arguments from mug in the translation-invariant case. Namely, 
implies 

Q(p ) 2 < 'P+(p)Q(p)'P+(p) - V°_(p)Q(p)V°_(p) 

and 

0 < tr C 4(|X>°(p)|Q(p) 2 ) < trc4 (T>°(p)Q(p)) 
for any p £ B( 0, A). Now, by Kato’s inequality, 


f \QW \ 2 

/R3 M 


dx<^ J tr C 4(|p| Q(p) 2 )dp < | J tr C 4(|r>°(p)| Q{p) 2 )dp 


for 


\P\ < \J\P? + ™o < P°(p)\ 

by ED: and therefore, when 0 < a < 4 / 7r, 


T(f) > T(/) + p(l-™/4) 


-8(0,A) 


tr<c 4 [|27 >° Cp) | (/ — f) 2 (p)]dp 


- (Wa,/4)m 0 ,, -,,2 

- ( 27r )3 V J\\ L 2(B(0 ,A),54(C)) • 


This immediately implies that / is the unique global minimizer of T and also shows 
that any minimizing sequence (/„) C A\ necessarily satisfies f n —> / strongly in 
L 2 (£(0,A),<S 4 (C)). 


• Step 3: regularity ofV°(p). It remains to show that the self-consistent free Dirac 
operator V° is smooth in the Fourier domain. To this end, we notice that it satisfies 
the equation 


V°(p) = D°(p) 


V°(q) 


4?r 2 Jb( o,a) \V 0 (q)\ \p~q\ : 


7 dq 


or, on .B(0, A), 


V° = D° 


a V° 1 

a^W\ * FF 


Therefore 

|VH>° = |V|D 0 -||1 

which shows that T>° £ H 1 (B( 0, A)). Applying now V to this equation and using a 
easy boot-strap argument, one obtains that 


V° £ p] H m (B{0,A)) cr(B(0,A)). 

m> 1 

which ends the proof of Theorem □ □ 
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3.3 Proof of Theorem 01 

Our proof proceeds as follows: we first minimize on a restricted set of translation 
invariant operators and study the limit of the corresponding sequence of minimizers 
as L —> oo. Using then the properties of the solution / constructed in Theorem^ 
we show that, for each L large enough, the so-obtained local minimizer is indeed 
the unique global minimizer of ££. 

• Step 1: Definition of the restricted minimization procedure. Let us define the 
following minimization problem 1 : 


El :=inf{££( 7 ) | 7 = /(-iV), / e B k } 


(75) 


where B\ is the set of functions / : r A i—> 54(C) which take the special form 
(compare with B\ defined in (1661) 1 



Vp e r%, f(p) = a ■ A(p) + (3b{p) 

(76) 

where A\T\ 

1—+ R 3 and b : r A 1—> R are such that 



Vp g r£, b(p) < 0, 

(77) 

Vperi, \A{p)\ 2 + b{ P f <\/a. 

The energy of such a state reads 

(78) 

££(7) 

= tr(£>°7 )-^/y \^(x 1 y)\ 2 W L (x-y)dxdy 

1 JJ(Tl ) 2 

(79) 


Y tr c 4 D °(p)f(p)) - 9r9 n 3 ^ 3 / \f(x)\ 2 W L (x)dx 

piS 2(2?r) Jt - 

Y tr C 4 (D°{p)f(p)) - a<y2n \ Y tr C *(f(p)f(q))W L (p - q), 


peri 


p.^erj 


where we have used that, due to the special form of /, the density of 7 vanishes, 
Pi = i_3 EperP tr c 4 f(p) = 0. 

• Step 2: Upper bound on the energy per unit volume. Let us now show that 

e t 

limsup— <E r = min{T(/) | / G ^l A } • 

L—>00 


This is easily done by considering the minimizer / of T as a test function. For the 
sake of clarity, we denote by the restriction of / to and by Gl = gL^—i^) 
the associated operator acting on £))(, which is the multiplication operator in the 
Fourier domain by (/(fc)) feer j- Notice that 

9 l(x) = £ f(k)e ik - (80) 

fcerj 

and that gL = /| r i belongs to B\ for any L due to the properties of / stated in 
Theorem ^ 

1 Recall that we denote also by —zV the operator acting on i.e. the multiplication operator 

in the Fourier domain by {k) ker l- Therefore /(—zV) is just the multiplication operator by 

(f(k))kerfc > i n the F° ur i er domain. 


25 






Using |/(fc)| < 1/2 for any k £ 13(0, A) and (IHHIl . one sees that 


\\9l\\i 


< 


(27r) 3 / 2 |rj 

2L 3 


\\V9lI 


< 


A(27T) 3 / 2 |rj 
2 L 3 


( 81 ) 


Therefore, is bounded in TU 1 , 0 O (K 3 ) and converges to / uniformly on compact 

subsets of M 3 . Finally, we have pa L = 0 since trc 4 f(k) = 0 for any k £ 13(0, A) and 
therefore, by C9, 

-4^ = 77 55 tT c*(D°(k)f(k)) - Q / \g L {x)\ 2 W L (x)dx. 

fcer^ ^ J • /Tl - 


Arguing like in the proof of Lemma ^ one easily shows that 


lim 

L —>oo 


£°l(g l ) 

L 3 


= nf), 


hence 

e t 

lim sup -kr- < T{f) = E t . 

L —>oo L 3 


• Step 3 : Existence of a minimizer for E]~. In order to show the converse inequality 

e t 

lim inf —> U r , (82) 

L—>oo L 3 ~ 

we first prove the existence of a minimizer for m- 

Since is finite-dimensional, there exists a solution 7 ° = /°(—*V), /° £ B 
of the minimization problem C3- We now argue like in the proof of Theorem |T| 
(step 1 ) to show that satisfies the same equation as global minimizers of £/ on 
Q\. Due to the fact that B\ is convex, f° also minimizes 

T /£(/)= E ^(V° L (p)f(p)) (83) 

P 6 Ia 


in B\, where V° L is the mean-field operator associated with f®=a- A° L + /? 


V° L (p) 


f°L(9)W L (q-p) 

a - (p - - 55 A °l(q) w l{9 ~ P)\ 

\ « er A J 

+P (m 0 - - 55 b °L(QW L (q 



Notice that 


to 0 - 


E b l^)W L {p ~9)> m 0 


(84) 


which shows that Tff (p) is invertible for all p. The unique global minimizer of T/ L 
on B\ is —sgn(T>^)/2, which is easily seen to belong to B](. We therefore obtain, 


for k £ T 


A’ 


/£(*) 


p £(*0 

2\V° L (k)y 


(85) 


26 












Saying differently, y£ := /°(—iV) satisfies the equation 


_.o sgn(£> 0 ) 

7 L =-o-> 


:= D° - ay ° L (x,y)W L (x - y ) 


or, written in terms of the projector V £ := 7 ° + 7^/2, 


V° L = X {- =o;0)(©£). 


• Step 4: Lower bound on the energy per unit volume. We now prove © by 
studying the weak limit of y£ = /°(—iV) defined in the previous step. Using the 
same type of estimate as (I%T1) . we obtain that (/°) is bounded in W 1,00 (R 3 ) and it 
therefore converges to some function /, uniformly on compact subsets of R 3 , up to 
a subsequence. On the other hand, let us define the stepwise function z/y,[/£] by 

MfM = E /°( fc )Xfc(p). (86) 

Xk being the characteristic function of n*=i[^*5 ^ + 27 r/L[n 3 ?( 0 , A). Then z/j,[/£] is 
bounded in L°°(B(0, A), 54 (C)) and thus J'lI/l] ~" h G A\, up to a subsequence, 
for instance for the weak topology of L 2 (B( 0, A), <S 4 (<C)). To see that h = /, we can 
for instance notice that 


/£(*) -£[/3(s) = ( 2^)- 3/2 E [ Xfc(p) - eip ' x ) + O ( 7 ) , 

where our 0(1/L) is due to the fceri which are close to the boundary of 73(0, A). 
Hence, for any compact subset K of R 3 , 

< C(27t)- 3/2 E l/°( fc )l [ Xk(p)\(k - p) ■ x\dp + O (j 

fcerj 7 V 


fl-Mfl] 


L°°(K) 


< %|B(0,A)|+O( -f 


for some constants C and C 7 , which shows that ^/-,[/})] —> / on compact subsets of 
R 3 and implies f = h G Aa- 
The energy of 7 £ reads 


^ = 73 E tr c^( 7 ?°(fc)/°(fc)) ^ 7W73 [ \fl{ x )? w L(x)dx 

L L ke r£ ^ ' Jtl 

+ |/xL 2 (p 7 J 2 . (87) 

We have 

lim f \f/ j (x)\ 2 WL(x)dx = f dx 

L^oo J Tl J R 3 |&| 

due to the fact that (/°)l converges locally to / and that 


7 i|lL»(T L ) 



( 2 ^) 3 |rj| 

4 L 3 
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is uniformly bounded in L. On the other hand, using the notation (PH) , 


j ^3 tr C *(D°(k)fl(k)) = ( 2 tt) 3 

ker A 


<B{ 0 ,A) 


tr C 4 (iy L [D 0 ]{p)u L [fl](p))dp + O 



Since vl[Il\ / weakly and vl[D°] —> D° strongly in L 2 (B( 0, A)), we deduce that 

r lim 73 H tr C 4 (-D 0 (^)/E(fc)) = (2 tt) _3 [ tr C 4 {D°(p)f(p))dp 
L ^°° L JB( 0 ,A) 

and 

tr C 4 (D°(k)fl(k)) - 2 ^ 3 J ^ |/£(a;)| 2 W L (a;)dx j =T(/). 


fcer4 


Recall that p, > 0, which implies due to E3 


f t 

liminf —> T(f) > E t 
L—>oo L 3 ~ K ’ ~ 


where we have used that / £ .4 a- This of course shows that 


pT 

lim = E 1 . 
L —>oo L 6 


and, as a matter of fact, T(/) = A 7 ". By Theorem^ / is the unique minimizer of 
T on .4a, and thus / = /. 

• Step 5: Study o/ i/ie convergence of 7 ° = /°(—zV). We know that —*• / 

weakly in L 2 (B( 0, A)) and that /£ —> / uniformly on compact subsets of R 3 . Notice 
that ^l[/ 3 ] £ -4a for any A. Since we know from the proof of Theoremnjthat T is 
continuous on .4a for the weak topology of L 2 (B( 0, A)), we obtain 

lim %[/“])=T(/), 


which means that (^l[/ 3 ])l is a minimizing sequence of T. By Theorem^we infer 
^lUl} f strongly in A 2 (£?(0, A)). Recall now that 

flip) = a ' A °l(p ) + Pb° L (p), f(p ) = « • u P fi{\p\) + Pfo(\p\) 

which easily implies that, more precisely, vl[A° l ] —> / 0 (| • |) and — * w-/i(| • |) 

strongly in L 2 (B( 0, A)). 

Using the equation PH fulfilled by /£, we are now able to prove a better con¬ 
vergence: 

Lemma 2. One has 

j = 0, (88) 

= 0, (89) 

j =0. (90) 

Hence, for L large enough, 

Vp£li, \V° l (p)\ 2 >\p\ 2 + ^. (91) 


lim \\V 0 L -V° l 

L —►oo 


ISoo(Tl) l 


„o „ Oil 


lim \\Vt - V°\ 


600 ( Tl ) l 


lim I 

sup 

\V°l(p) ~ 

-V°(p) 

L—>o 0 ' 

^peri 



lim 

sup 

\f°(p) - 

- m 

L—>00 

\P eI A 


j 

lim I 

sup 

I K(p) - 

-v\p)\ 

L —»oo < 

^per^ 
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Proof. Since easily implies El, ® and (El) , we only have to prove that 


sup 

pe n 


(27r) 5 / 2 7T 

L 3 


E A° L (p-q)W L (q ) 


w p _ g /i(|p- q\) 


'B(O.A) 


kl 2 


dg 


and that 


(92) 


sup 

P^A 


(27r) 5 / 2 7T 

L 3 


E b °L(p-Q)W L (q) 

1^ L a 


'B(O.A) 


/o(b- g|) 

kl 2 


dg 


(93) 


as L —> oo. We only treat El) : El is obtained by the same arguments. This is 
easily done by noting that for instance 


(2tt) 5 / 2 7 


L 3 


E b L(P~ l) w L{q) - f 

Je 


?er- 
l<?l<p 


/o(b- g|) j 

I# * 


< 


B( 0 ,rj) 


2 n 2 p f 2 -k\ 

— + u) 


\ 


E 


V 


<?er£ 

9#0. |g|<7j 


1 

kf 


'B(0 


dq 


/ 


due to 16° (p)| < 1/2 and \fo(p)\ < 1 / 2 , and an easy estimate for the rest of the 
form 


(2tt) 5 / 2 7 


L 3 


E b °L(P-q)W L (q)- I 

Je 


\q\>v 


/o(b-g ) j 
- m 2 -«g 

B(0,A)\B(O.J7) 191 


< 


VL[b Q L \{p~q)v L [ 1/| • | 2 ](g) - ^ ) dg 


0(1/L) 


I B(0,A)\B(0,v) \ kP 

— ^ IM*>£] - /o||l 2( B ( 0jA )) + C II[1/1 • | ] - 1/1 • I Hi,2(b(o,a)\b(o,?7)) + 0(1/L) 

which converges to 0 as L —> oo, due to the strong convergence of to /o in 

L 2 (B( 0,A)). □ 


• Step 6: 7 ° = /£(—*V) zs f/ie unique global minimizer of £° on Qfor L large 
enough. Using El , we now show that 7 ° = V 3 — lfc/2 is the unique global 
minimizer of £°. Indeed, using again ideas from □ G 2 we write for any 7 e Gk: 

£°l( 7) = SK'rD+tviViQ)-^ ff | Q(x,y)\ 2 W L (x-y)dxdy 

1 J J(Tl) 2 

+ ^D l (pq, pq) 

where we have used that, by construction, p^o(x) = 0 for any x £ R 3 , and where 
( 5=7 — 7 °. Notice that Q satisfies the inequality 

-V° L <Q< 1-V° L 

which now implies as in El EH 

0 <tr(|P°|Q 2 )<tr( 2 ?°Q). 
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By Lemma ^ we can estimate 


a 

2 


ff \Q( X , y)\ 2 W L (x — y)dx dy < ^-Cfc(m 0 /2) tr(J-A + ml/AQ 2 ) 

z v 


< 2 C A(mo/ 2 )tr(|P°|Q 2 ) 


due to ED- and therefore 


tr (D° L Q)-^ff \Q(x, y)\ 2 W L (x — y)d. 


- y)dx dy > 1 - 


aCf (toq/ 2 ) 


On the other hand, Dl(pq, Pq) > 0 and therefore 

> £l(il) + (i - aCk( ™° /2) 

For L large enough one has 




tr(|2?£|Q 2 ) 


(94) 


1 _ aCjj{jno/2) > Q 


since 0 < a < A/it and limA.,^ C^(mo/2) = Ca (mo/2) < 7r/2 by Lemma Q 
Therefore, by Lemma [3 and 1941) . 


mi) > mil) 


mo (l - ctCa (too/2)/2) 


7l - 7 




which implies that y£ is the unique minimizer of ££ on Q^ and ends the proof of 
Theorem H J □ 


3.4 Proof of Theorem [5| 

For the sake of simplicity, we only treat the case where 

H| = ^(P°k+ ^ W L { Xi - Xj ) 

i> 1 1 <i<j 

is chosen instead of H£, the arguments being easy to extend to H£. 

The energy of a Hartree-Fock state of density matrix O (a self-adjoint operator 
acting on such that 0 < 0 < /£) is now simply f£(0) instead of (BUI) (recall 
that ££ is defined in ED)- Therefore, restricting this energy to translation-invariant 
operators, the minimization problem to be solved for a fixed L is 

E l :=min{££(0) | 0* = 0, 0 < 0 < jf, | 0 = f(-*V)} . (95) 

Since 53 a i s finite dimensional, there exists a minimizer Ql(x, y) = (2Tr)~ 3 ^ 2 ^ L (x—y) 
of this energy, where 


&(*) = m/L) 3/2 0 L (k)e k (x). 


Computing the energy of 0^, we obtain 

m&L)= tvc4D°(k)e L (k))-L 3 -^ [ |a(a:)| 2 W L {x)dx+^pL\p &L ) 2 . 

^ 2{27T) 2 

(96) 
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Notice that 


E tr C 4 (D°(k)Q L (k)) 


fee ri 


< V A2 + m o l r A< 


since 0 < ©i(fc) < /c 4 for any k, and that 

[ \f L (x)\ 2 W L (x)dx < Ck(mo)(\D°\Z L ,t L ) 

JTl 


= E \D°(k)\tv c <(Q L (kf) 

fee 


< 4(27t) 3 C'a(too) y^A 2 + m, 


2\£a[ 

L 3 


by Lemma 0 This implies 


L 


-3 


E tT C 4D°(k)Q L (k))-L 3 -^— / 

2 ( 27r ) JT l 


feer^ 


< M 


for some uniform constant M. Hence, using £°(©z,) < 0, we obtain from lIHHIl 

\lxL\p &L f<M 

and limL-xx, p@ L = 0. To end the proof of Theorem 0 we notice that 

( 2 t r ) 3 / 2 

fee rj 

< y: tr C . e L (t) 


ia I 


L °°( Cl ) 


ll£felL“(R 3 ) — j^3 E 


feerj 


< C(2n) 3 / 2 pe L 

where we have used that 0 l > 0 and that \M\ = M) 1 / 2 < Ctrc 4 (|M|) for 

some constant C. This proves (15^1) . □ 

Remark 4. Passing to the weak limit in the energy, it is also easy to prove that 
limL^oo El/L 3 = 0. 


3.5 Proof of Theorem EH 

• Step 1: Construction of a minimizing sequence of finite rank operators for £g DF . 
In order to prove an upper bound for the energy difference El(p) — El( 0), we need 
to use a trial state. To simplify matter, we shall use a finite rank operator. 

Lemma 3. Assume that 0 < a < 4/7r, A > 0, mo > 0 and that n £ C. Then 
there exists a sequence (Qk) C Qa of finite rank operators such that Qk +V°_ is a 
projector for any k, and 


lim <?b DF (Q/c) — min£g DF . 

k— «x> Qa 

Proof. We know from Theorem [T] that £g DF possesses a minimizer Q = V- — £ 

Qa such that V- is a solution of the self-consistent equation 

P- = X( —oo;0) iP° + PQ ~ Rq) 
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where 


/ X 1 O Q(*»?/) 

VQ ■= u{PQ - n) *—y R Q- =a ^—^y 

We know from 22,H.2J2 that Vi pq £ L 2 (R 3 ), yjg £ i 6 (R 3 ) and i?g £ 6 2 C£j A )- Now, 
let be (tpk) and (Rk) two sequences such that 

1 . ipk £ i 1 (M 3 ) niJ 1 (M 3 ) for any k, Vifik —► Vy>g in L 2 and ^ —> yjg in L 6 (R 3 ) 
as k —> oo; 

2 . i?fc £ Si(^a) for any k and Rk —► in 0 2 (7} A ) as fc —> oo. 

We define P k := X(-oo; 0 ) (2?° + <Pk - Rk) and Q' fe := P k - V°_. For the sake of 
simplicity, we assume that 0 is not in the spectrum of V° + <pq — Rq , the following 
arguments being easily adapted to the other case. Then 0 ^ a(V° + ip k ~ Rk) for k 
large enough and we can use Cauchy’s formula like in [2171 [25) to obtain 

1 1 

T>° + <pk — Rk + if] + if] 

l r°° i i 

= n / ^oT- (y>fc ~ 24) ^ . dr] (97) 

2-7T J_ OQ 2? u + ip k - Rk + IT] V U + IT] 

which implies Q' k £ ©i(5j A ) for 




1 


V° + ip k - R k + it] 


(ipk - Rk) 


1 


V° + ir] 


6i(Sa) 


< 


\/e 2 + if sj m,Q + if 


II Pk ^fcllsi (jjy, 


since | 2 ?° + ipk — Rk\ > e > 0 and thanks to lT7Tl) . Notice that ipk £ Si(^a) since, 
denoting by 1 a the characteristic function of the ball B( 0, A), 


I^ fc Isi(J5 A ) “ 


< 


|l A (-iV)^fcl A (-*V)|| Si(j5A) 
l A (-iV)|v5 fe | 1/2 2 ^ 


^ c-IIIaII^ II^ILi 


by the Kato-Seiler-Simon inequality, see m and )51 Theorem 4.1], or the appendix 
of 121 ) . 

We now show that lim*,-^ £g DF {Q'k) = £g DF ((3). D ue to th e results of Klaus- 
Scharf jH3J or the estimates of )25], we already know that Q' k —> Q in 6 2 (#a)- Since 
£bdf i s strongly continuous ESI for the norm 


II <3II == (llQ| 2 6a ( flA ) + IIP 0 lc) 1/2 . 


(98) 


it therefore remains to prove that pQ' k —> pq in C. Notice that due to the ultraviolet 
cut-off, pQ' k —> pq strongly in L 2 (see, e.g., 23 Formula (9)]). Expanding llbTI) . we 
obtain 


^=-^e/ ^(n 


j -1 ■ 


1 V° + in 

i =1 


— (74 - <p fc ) 


V° + irj 


*/>( ]M 


z?7 


~{Rk Pk) 


1 


2 ?° + V^fc - 74 + if] 


(99) 
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As proved in [221 Section 4.3.3], the first sum of the r.h.s. of 1991) is continuous for 
the norm (19811 when \7ipk —► V<^>q in L 2 (R 3 * ), ip k —> Pq in T 6 (R 3 ) and Rk —■> Rq 
in e 2 (i3A) (notice that the proof of [231 has to be adapted to the case where D 0 
is replaced by T>°, which is an easy task). It therefore remains to show that the 
density associated with the last term of (PI) is also continuous for the C norm as 
k —> oo. The only non trivial term is 

dl1 (n V° + V° + ipk - R k + irj 

which indeed is continuous for the Si(f)A) topology by estimates of the type 


n 

U =1 

< c 


1 


1 


T>° +irj^ k I T>° + ifik ~ Rk + irj 

6 i 


6i(6a) 


1 


< c 


V° (p) + irj 

1 


hkf L e 


D°(p ) + irj 


lb 

6 


L 6 


\J e 2 + V 2 


\Je 2 + rf 


= C(m 2 0 + V 2 )- 3 / 2 (e 2 +r, 2 )- 1 / 2 (1 + p 2 )~ 1/2 L Jp k \\ 


6 

L e ' 


This means that the associated density (p k ) converges in L 1 . Therefore the Fourier 
transform (p)T) converges in L°° and has a compact support by assumption. Hence 
(p k ) converges in C. 

As a conclusion, Q' k —s > Q as k —> oo for the norm ||Q|| := (|Q|| e2 (^ A ) + IIpqIIcT 2 - 
It is proved in [231 Section 5.1] that £g DF is strongly continuous for this norm, 
which implies now lim k ^-oo ^BT>F(Q k ) = £g DF (Q)- Hence, we have constructed 
a minimizing sequence of trace-class operators in Q\. As a last step, we now 
approximate each Q' k by a sequence of finite rank operators (since k is fixed, we 
shall denote for simplicity Q' = P' — P^_ instead of Q' k ). 

This will be done by using the following decomposition 2 proved in m 


N 


P'=>\fn)(fn\ + 


\Ui \iVi) (ill ~\~ XiVi\ 


n= 1 


i= 1 


1 + A f 


( 100 ) 


1 -P' = 


M 

E 

m— 1 


1 9m) (gn 


OO 

E 

2=1 


|Vi - \iUi){vi - XiUi 

TTa 2 


( 101 ) 


where U (vi)i>i is an orthonormal basis of P°Sja, (ffi)f=i U (uj)i>i is an 

orthonormal basis of P°Sj a, and JT^A 2 < oo. Using m and tnm , one sees 
that 


N M . 2 

Q' = E “ E \9m){gm\ + E ! + \2 (l^iX^I ~ l«*X“t|) 

n—1 m =1 z>l * 

+ E 1 + + KX^l) ( 102 ) 

i> 1 * 


2 The reader should compare this decomposition with the classical formula which gives the BDF 

state associated with the projector P ^ in the Fock space built with see, e.g., EH Formula 

(10.96)] where A = Aj|«i)(vi|. 
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which easily implies that a(Q') fl ( 0 ; 1 ) = (Ji>i { i 




. Recall that Q' £ 


®iC£)a) which implies (A,;), £ O(R). Let us now define 

N K 


r> i c \ / f I , \ Ui + ^ i Vi ) ( U i + I , I \ / I 

P K ~ 2^\M(fn\+2_^ -j—Ti- + 2^ \ U U\ U i I' 


i =1 


i=if+l 


The operator 


AT M K \2 

Qk = PK -V°_ = E \fn)(fn\ - E \9m){g m \+Y x + \2 (l^X^I ~ KXfol) 
n=l m=l 2—1 z 

+ X! 1 + A2 (K)(^l + ki)(^l) ( 103 ) 

2=1 * 

is a finite rank operator in Qa such that Qk —> Q' for the Si(f)A) topology as 
K — > oo. Therefore lim/c^oo £b DF (Qa") = ^bdf (Q')- This ends the proof of Lemma 

El □ 


• Step 2 : Upper bound on the energy difference El(<p) — E^{ 0). Let be e > 0 and 
Q = P — Vf £ Qa a finite rank operator of the form for some K 1 such that 

i ^bdf( ( 5) — niinfgoF + e - 


We now define, for i = 1, K , 


«?(*) = E ^(k)e ik '* 

kerfc 

and (fn)n= i) (Sm)m=b ( v i)iLi by similar formulas. We have, using the projector 
■p£ defined in Theorem El 


(p£nf,P°^ L ) i2(Ti) = ('=pj E (^°(fcKi(fc),%(fc)) c4 

' ' fcerj 

and so 

lim {V^uf,V^uf) = [ (V°_(k)u i (k),u j (k)) ci dk = 5i j . (104) 
On the other hand, we also have 


lim 


- V 


0 „,£| 


Wl 2 (T l ) 


= (v: 


-f 1 


'«A 


= 0. 


(105) 


Using the same type of behavior for (/7)2=i> (Sm)™=l an( i ( v i)iLi and the Gram- 
Schmidt orthonormalization procedure, we may therefore find an orthonormal sys¬ 
tem (/^)2= i U (57)^i of P^a and an orthonormal system (<?£)£f =1 U (uf')^ 1 of 
(1 — Va)^>a suc h that 


lim II uf — uf 

L—>-oo 


£ 2 ( Ti) °> 


lim flvf-vf 

ly—KX> 


A 2 (Ti,) 0 


(106) 


for all i = 1 ,..., iv, and 


lim 


f - r 

«/ 72 J n 


A 2 (Ti) 


= 0 , 


lim 


T?m||z,2(T t ) 0- 


(107) 
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for all n = 1, N and m = 1, M. 
We now define our trial state by 


N 


M 


K 


Ql = Yl \fn)(fn \ - IflmXffml + 1 + \A \ d i)^i\ ~ 


K 


+ Y1 l + ‘x2 + \Vi){Ui\) ( 108 ) 

2=1 1 


and 7 l := 7 ° + Ql, where we recall that 7 ° = — 1^/2 is the unique translation- 

invariant minimizer of defined in Theorem 0J Notice that by construction, Ql 
satisfies 

-V° l <Q l < 1 - Vl 

and therefore 7 l & G\- Let us now compute 


££(%)-E L {0) = SWl)-^l) 

0 


= Wl+Ql)-S till) 

tr {V° l Ql)- [ <Pl(x)pq l (x)(Ix+?-D l (pq l ,Pq l ) 

* 


(T l ) 2 


\Ql(x,v)\ 2 Wl(x - y)dxdy , 


(109) 


where we have used that p 1 0 = 0. Passing to the limit in (110911 and using ll(Hil) and 
(fTTTTIi . one easily obtains 


lim (f£(7 L ) - E l ( 0)) = £bdf(Q) < min^DF + e 

L—>oo Qa 


and therefore 


lim sup(-El(</?) - E l ( 0)) < min£g DF . 

L-> 00 Ca 


• Step 3: Lower bound on the energy difference El{<p) — El(0). Since is finite- 
dimensional, there exists for any L a minimizer j' L of £f on Q\. Using classical 
arguments already used in the Hartree-Fock theory HE El EH (see also |26l Lemma 
2 ]) and the positivity of Wl (except on a set of measure zero), one easily shows that 

+ La/2 = Vl + A|yj)(</?| 


where Vl is an orthogonal projector, A € [0; 1] and ip £ ker(V 7 > ) with 
X>y L = D° + a(py L - n L ) *W L - al L (x,y)W L (x - y) 

Since then £1(Vl — ijf) = £lI'l)i 7 l '■= Vl — I \/2 is also a global minimizer of 

£t on Ql 

Let us now define Ql := 7 l — 7 ° = Vl ~ V°. Like in 11091) . we have 


E L (ip)-E L (0) 


£IIl)-£HI) 
tr (V° l Q l ) - [ 



a 

Pl{x)p Ql (x)dx + -D l (pq l ,Pq l ) 
\QL{x,y)\ 2 W L {x - y)dx dy. 


(110) 
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Therefore, arguing like in [3 ESI ESI and using lEHT) and (I7TT>1) . we infer 


E L (<p)-E L ( 0) > (l^aCk(m 0 /2)/2)tr(\V° L \(Q+ + -QZ~)) 

a a 

+ 2 Dl (PQl ~ n L, PQ l ~ n L ) - -D L (n L , n L ) 

> {l-aCl{m 0 /2)/2)tr(p 0 L \Ql) 

CM. a 

+^ d l{pq l - n L , pq l - n L ) - - D L (n L , n L ) (111) 


where Q^ + := (1 - V®)Ql( 1 — P°) > 0 and Q2~ ■= 'PlQl'Pl < 0- Since 
lim^^oo C^(m 0 /2) = Ca(too/ 2 ) < it/2 and 0 < a < 4/7r by assumption, then 
1 — aC^(m 0 /2)/2 > 0 for L sufficiently large. Using tr (|X> 3 |Q 2 ) > m 0 /2tr (Q\) 
due to ct . we therefore deduce from cm that 

1. QL(x,y)tc L { x )tc L (y) is bounded in L 2 ((R 3 ) 2 ) and therefore Pq l (x)1c l ( x ) i s 
bounded in L 2 (R 3 ) ; 

2 - P|x>o|i/ 2 q++|-do p/ 2 (a:)lei, {x) andp_ |:D o |1/2 Q £ - |2 ,o|i/ 2 (a;)lc i (a:) are non-negative 
functions, uniformly bounded in L 1 (K. 3 ). 


Thanks to the ultraviolet cut-off, we may therefore assume that 

1 . QL(x,y)tc L ( x )tc L (y) — 1 Q(x,y) weakly in L 2 ((R 3 ) 2 ) and uniformly on com¬ 
pact subsets of R 6 ; 


2. pq l (x)tc L ( x ) ~ 1 Pq( x ) weakly in L 2 (R 3 ) and uniformly on compact subsets 
of R 3 ; 


3 - P|r.o|i/ 2 Q ++|i 50 |i/ 2 (x)l CL (x) -> p+(x) £ T 1 (R 3 ) and p_ |x ,o |i /2QZ -| C o |i /2 (a) 
1 c L {x ) —> p~(x) £ T 1 (R 3 ) uniformly on compact subsets of R 3 . 


Let us first show that Q £ Q a. Due to ro . one easily obtains, passing to the 
weak limit, 


-V°_ <Q<V% 


Notice also that \V° L \ 1 ^ 2 Q L | x / 2 (a;, y) = \D° L \ 1 / 2 'P^Ql'P { 1\'D {) l \ 1 / 2 {x, y) converges 

uniformly on compact subsets of R 6 to | D°\ 1 / 2 'P ( ?_Q'P ( ?_\'D 0 \ 1 / 2 (x,y), by (IHHll and 
Therefore, one obtains 


P + — P|D°| 1 /2Q++|-DO|l/2 £ L^R 3 ), p — p_| X )0|l/2Q-|D0|l/2 £ L 1 (R 3 ) 

where, this time, = V^_QV^_ and Q ++ = V+QV+, which shows that Q is 

_ v o _____ 

P^-trace class, Q £ S 1 “(1)a). Finally, 11111) shows that D l (pq l — ul,Pq l — til) 
is bounded. Passing to the limit, we obtain pq — n £ C which implies pq £ C and 
as a conclusion Q £ Q\. Hence, it remains to show that 


liminf (. E L {ip) - E L (0 )) > £^ DF (Q)- 

L —KX) 


This will imply the desired bound 

liminf ( E L (ip ) - E L { 0)) > min££ DF 
L—>oo Qa 


since Q £ Q\. 

By Fatou’s lemma, we have 

liminf D l {pq l - n L , Pq l - n L ) > D(pq -n,pq- n) 

L —>oo 
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and since obviously limL^oo Dl(ul, til) = D(n,n) due to the approximation of 
integrals by Riemann sums for the continuous function n, it only remains to prove 
that 


liminf (tr (\T>° L \(Q+ + - Q L )) - ^ [f \Q L (x,y)\' 2 W L (x - y)dx dy 
L ^°° \ i J J(J L ) 2 

> tr(|2?°|(Q ++ - Q~~)) ~ 77 [[ ^ X ' V ^ dxdy. ( 112 ) 

2 e \x - y | 

This is shown by following exactly the method the authors of m used for the proof 
of their Theorem 1 (step 1), which we briefly outline now. 

The idea of fSB] is to use space cut-off functions ijr and defined by iir(x) = 
y(\x\/R) and £r(x) = £(|x|/.R) where 77 , £ £ C°°([0; oo); [0; 1]) are such that r ? 2 +£ 2 = 
1, 77 (f) = 1 if t £ [0; 1] and 77 (f) = 0 if f > 2. Let us denote by the periodized 
cut-off function which is defined for L large enough by the same formula as tjr on 

C L , and by = yjl- (y^) 2 . 

Then, in order to reproduce the first step of the proof of |2S1, Theorem 1], the 
following analogue to (23 Lemma 1] is needed: 

Lemma 4. We have 


lira limsup || [|X>i|,^] 

00 l ^ oo 


SocCSj 


= 0 . 


Proof. We have 

II W^Mr] 


<2 II \V° L \-\V°\ 


[l^°l ,&] 


(113) 


le«(*J;) - " II '^1 “ W ' Il6„,(fli) ^ II l>«J IICooWJ:) ' 

The first term of the right hand side of m tends to 0 as L —> 00 , due to Lemma 
13 By the regularity of T>°(p ) proved in Theorem 0 there exists a constant C such 
that 

\\V 0 (p)\-\V°(q)\\ <C\p-q\ 

when p, q £ Tj(. This enables us to argue similarly to the proof of El, Lemma 1] 
and obtain a bound of the form 




|| [|Z>°I ,&} 

and therefore 

limsup I [|T>°Uh] || s (fli) 

L—>00 

which ends the proof of Lemma 01 
One has 


C v- 
S 75 £ 


p£r(p) 


P £(2ttZ. 3 )/L 


<c \r^ R (r)\dr = 0(l/R), 
Jr 3 


□ 


Ml v l\(Ql + -Ql )) = ^\V° l \(Q++ -Q l )t}r) 

+ tr(| -D° l \^(Q++ - Ql~)&) + tr([& \V° L \}(Q++ - QZ~)&), 

On the other hand, 


[[ \Ql(z.v)\ 2 Wl{x — y)dxdy = [[ 

J J(c L ) 2 JJci \x-y\ 

+ JJ 2 (£r) 2 ( x )\Ql(x, y)\ 2 W L (x - y)dxdy + 0(1/L) + 0(1/R) 
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where we have used 116411 . Now, multiplying — V® < Ql < 1 — V £ by on both 
sides and using Kato’s inequality m for L large enough, we easily obtain, following 
[2II], that 

trClX’il Zr(Ql + ~ Ql~)£,r) ~ j [({£,R) 2 {x)\Q L {x,y)\ 2 W L {x-y)dxdy > 0. 

^ j Jc 2 , 


Hence, 


t r (I^Ll((3i + ^ Ql )) - \ JJ^\Q L {x,y)\ 2 W L (x-y)dxdy 

V 2 R ( x )vl R (y)\QL(x,y)\ 2 


>^(Vr\VlKQl -Ql ) 7 7fl)-^ //c2 


\x-y I 


dy 


+ tr([&, |^|](QJ + - QZ")&) + 0(l/i) + 0(l/ii). (114) 


The estimate m is then easily obtained by passing to the limit first as L —> oo 
and then as R —> oo, and using Lemma ^ 

As a conclusion, we have shown that 


lim ( E L (ip) - E l ( 0)) = minfg, 

L— >oo Qa 


Due to our estimates, we also conclude that the weak limit Q of Ql satisfies 
^bdf (Q) = m i n QA^BDF- Therefore, Q is a minimizer of 5 bdf- This ends the 
proof of Theorem [SJ □ 
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